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We report on a search for gravitational waves from coalescing compact binaries using LIGO and 
Virgo observations between July 7, 2009 and October 20, 2010. We searched for signals from binaries 
with total mass between 2 and 25 Mq; this includes binary neutron stars, binary black holes, and 
binaries consisting of a black hole and neutron star. The detectors were sensitive to systems up to 
40Mpc distant for binary neutron stars, and further for higher mass systems. No gravitational-wave 
signals were detected. We report upper limits on the rate of compact binary coalescence as a function 
of total mass, including the results from previous LIGO and Virgo observations. The cumulative 
90%-confidence rate upper limits of the binary coalescence of binary neutron star, neutron star- 
black hole and binary black hole systems are 1.3 x 10"", 3.1 x 10~® and 6.4 x 10~® Mpc~^yr~^, 
respectively. These upper limits are up to a factor 1.4 lower than previously derived limits. We also 
report on results from a blind injection challenge. 



PACS numbers: 95.85.Sz, 04.80.Nn, 07.05.Kf, 97.60.Jd, 97.60.Lf, 97.80.-d 
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I. INTRODUCTION 

During 2009 and 2010, both the Laser Interferometer 
Gravitational-wave Observatory (LIGO) [1] and Virgo [2] 
gravitational-wave detectors undertook science runs with 
better sensitivity across a broader range of frequencies 
than previously achieved. Among the most promising 
sources of gravitational waves for these detectors are 
compact stellar mass binaries as they spiral in toward 
each other and merge. For such systems, which include 
binary neutron stars (BNS), binary black holes (BBH), 
and neutron star-black hole binaries (NSBH), the late 
stages of inspiral and merger occur in the most sensitive 
band (between 40 and 1000 Hz) of the LIGO and Virgo 
detectors. In this paper, we report on a search for grav- 
itational waves from binary systems with a maximum 
total mass of 25 Mq , and a minimum component mass of 
IMq. 

A hardware injection was performed during the data 
collection without the knowledge of the data analysis 
teams as part of a "blind injection challenge" [5]. This 
challenge was intended to test the data analysis pro- 
cedures and processes for evaluating candidate events. 
The injection was performed by coherently actuating the 
mirrors on the LIGO and Virgo detectors to mimic a 
gravitational-wave signal. Prior to its unveiling as an in- 
jection ( "unblinding" ) , the event was determined to be 
a candidate gravitational wave: it was found to have a 
false alarm rate of less than 1 in 7, 000 years and no evi- 
dence for an instrumental or environmental origin could 
be found. After the analysis of the event was finished it 
was revealed to be a blind injection and removed from 
the data. 

With the blind injection removed there were no gravi- 
tational waves observed above the noise background. As 
a result we place upper limits on rates of compact bi- 
nary coalescence (CBC), using upper limits from previ- 
ous gravitational-wave searches [4] as prior information. 
The upper limits presented here are up to a factor 1.4 
lower than previously derived limits but still two to three 
orders of magnitude above expected CBG rates [5]. 

The paper is laid out as follows. In Section [Ilj we pro- 
vide a brief description of the detectors and their sensi- 
tivities during LIGO's sixth science run (S6) and Virgo's 
second and third science runs. In Section [lll| we present 
a brief overview of the analysis methods used in perform- 
ing the search. In Section |IV| we discuss the recovery of 
the blind injection. In Section|V]we present the results of 
the search with the blind injection removed. In Section 
|VI| we give the upper limits obtained from the search and 
close with a brief discussion in Section IVlIl 



II. DETECTORS 

The LIGO observatory comprises two sites, one in Han- 
ford, WA and the second in Livingston, LA. The data 
used in this search were taken during S6, which took 




Frequency (Hz) 



FIG. 1: Typical detector strain noise spectral density for the 
LIGO S6 and Virgo VSR2/3 runs. From lowest to highest at 
10^ Hz, the curves are for the HI, LI and VI detectors. 



place between 7 July 2009 and 20 October 2010. Dur- 
ing S6 each of these sites operated a single 4km laser 
interferometer, denoted as HI and LI respectively. The 
2km H2 instrument at the Hanford site which operated 
in earlier science runs was not operational in S6. Fol- 
lowing LIGO's fifth science run (S5) [1 , several hard- 
ware changes were made to the LIGO detectors so that 
prototypes of advanced LIGO technology could be in- 
stalled and tested [51 [7] . This included the installation 
of a higher power laser, and the implementation of a DC 
readout system that included a new output mode cleaner 
on an advanced LIGO seismic isolation table 8j . In addi- 
tion, the hydraulic seismic isolation system was improved 
by fine tuning its feed-forward path. 

The Virgo detector (denoted VI) is a single, 3km laser 
interferometer located in Cascina, Italy. The data used in 
this search were taken from both Virgo's second science 
run (VSR2), which ran from 7 July 2009 to 8 January 
2010, and its third science run (VSR3), which ran from 11 
August 2010 to 20 October 2010. In the period between 
the first Virgo science run (VSRl) and VSR2, several 
enhancements were made to the Virgo detector. Specifi- 
cally, a more powerful laser was installed in Virgo, along 
with a thermal compensation system and improved scat- 
tered light mitigation. During early 2010, monolithic sus- 
pensions were installed, which involved replacing Virgo's 
test masses with new mirrors hung from fused-silica fibers 
[S]. VSR3 followed this upgrade. 

The sensitivity of the detectors during the S6, VSR2 
and VSR3 runs is shown in Figure [ij The correspond- 
ing sensitivity to binary coalescence signals is shown in 
Figure |2] This figure shows the distance at which an 
optimally oriented and located binary would produce a 
signal-to-noise ratio (SNR) of 8 in a given detector. The 
figure illustrates the improvement in sensitivity for the 
LIGO detectors between S5 and S6 and for Virgo be- 
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Although the template waveforms in this search neglect 
the spin of the binary components, the search is still capa- 
ble of detecting binaries whose waveforms are modulated 
by the effect of spin PTI . 

We require candidate signals to have a matched filter 
SNR greater than 5.5 in at least two detectors, and to 
have consistent values of template masses and time of 
arrival (allowing for travel-time difference) across the de- 
tectors where this threshold is exceeded [TSj. We use a 
chi-squared test I19j to suppress non-Gaussian noise tran- 
sients, which have a high SNR but whose time-frequency 
evolution is inconsistent with the template waveform. If 
the reduced chi-squared of a signal, Xrj is greater than 
unity, we re-weight the SNR p in order to suppress the 
significance of false signals, obtaining a re- weighted SNR 
statistic"'^ 



FIG. 2: Inspiral horizon distance versus the total mass 
of equal-mass binaries from S5/VSR1 (gray lines) and 
S6/VSR2/VSR3 (colored lines). The horizon distance is the 
distance at which an optimally located and oriented binary 
would produce an expected signal-to- noise ratio of 8. The 
figure shows the best sensitivity achieved by each detector 
during the runs. 



tween VSRl and VSR2. The reduction in the horizon 
distance of the Virgo detector in VSR3 is due to a mir- 
ror with an incorrect radius of curvature being installed 
during the conversion to monolithic suspension. 



III. BINARY COALESCENCE SEARCH 

To search for gravitational waves from compact binary 
coalescence [H [TOl E] , we use matched filtering to corre- 
late the detector's strain output with a theoretical model 
of the gravitational waveform [12]. Each detector's out- 
put is separately correlated against a bank [T31 of tem- 
plate waveforms generated at 3.5 post-Newtonian order 
in the frequency domain |14[ I15j . Templates were laid 
out across the mass range such that no more than 3% 
of the SNR was lost due to the discreteness of the bank. 
Only non-spinning waveforms with zero eccentricity and 
a component mass > I were generated, and the tem- 
plates were terminated prior to merger. In the early 
stages of the run, as in previous searches [H [101 [H] i the 
template bank included waveforms from binaries with a 
total mass M < 35 Mq . However, the search results indi- 
cated that the higher mass templates (M > 25 Mq) were 
more susceptible to non-stationary noise in the data. Fur- 
thermore, it is at these higher masses where the merger 
and ringdown phases of the signal come into the detec- 
tors' sensitive bands. Consequently, the upper mass limit 
of this search was reduced to 25 during the latter 
stages of the science run. Results of a search for higher 
mass binary black holes using non-spinning, full coa- 
lescence (inspiral-merger-ringdown) template waveforms, 
such as in iI6j , will be presented in a future publication. 



[(1 + (X^)^)/2]1/^ 
P 



for Xr > 1: 

for Xr < 1- 



(1) 



Our analysis reports the coalescence time and the 
quadrature sum, pc, of re- weighted SNRs for events 
coincident between the detectors. The statistic pc is 
then used to rank events by their significance above 
the expected background. To measure the background 
rate of coincident events in the search, we time-shift 
data from the detectors by an amount greater than the 
gravitational-wave travel time difference between detec- 
tor sites and re-analyze the data. Many independent 
time-shifts are performed to obtain a good estimate of the 
probability of accidental coincidence of noise transients 
at two or more sites. The analysis procedure described 
above is similar to the one used in previous searches of 
LIGO and Virgo data, such as [TD] and [IT]; it will be 
described in more detail in [20] . 

The background rates of coincident events were ini- 
tially estimated using 100 time-shifted analyses. These 
background rates vary depending on the binary's mass 
— via the waveform duration and frequency band — 
and also on the detectors involved in the coincidence (the 
event type). The relevant mass parameter is the binary's 
chirp mass, A4 = {m,im2)^^^{'mi + 7712)"^^^, where toi 
and TO2 are the component masses in the binary system. 
Thus, we sort coincident events into three bins by chirp 
mass, and by event type |10j . 

The requirement of a coincident signal between at least 
two sites restricts the times that can be analyzed to 
four distinct types of coincident time. Between July 
2009 and October 2010, a total of 0.56 yr of two-or- 
more-site coincident data was collected. This comprised 



Equation [T] is an improvement over the "effective SNR" used 
to rank events in [4l[l0l[TT]. Most notably: while effective SNR 
also re- weighted SNR using Xri It became larger than SNR when 
< 1- This made it susceptible to over weighting events that 
had statistical downward fluctuations in yj. 
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0.14 yr of HILIVI coincident data, 0.21 yr of HlLl data, 
0.13 yr of HI VI data, and 0.08 yr of LI VI data. Dur- 
ing HlLlVl coincident time there are four distinct event 
types: HILIVI, HlLl, HlVl, and LlVl. In S6/VSR2, 
all four event types were kept. In S6/VSR3, HlVl and 
LlVl events in triple-coincident time were discarded due 
to the heightened rate of transient noise artifacts in Virgo 
and its decreased sensitivity. 

For each candidate, a false alarm rate (FAR) is com- 
puted by comparing its pc value to background events 
in the same mass bin and coincident time and with the 
same event type. Candidates' FAR values are then com- 
pared to background events in all bins and event types, 
over the appropriate coincident time, to calculate a com- 
bined FAR. This is the detection statistic which is used 
to assess the significance of events over the entire analysis 
time. 

Due to the finite number of time-shifts performed, the 
smallest non-zero FAR that can be calculated is 1/Tbg, 
where Tbg is the total background time obtained by sum- 
ming the coincident live time in each time-shift. If an 
event was found to be louder than all background events 
within its analysis period, additional time-shifted analy- 
ses were performed to calculate a more precise FAR for 
the event. 

Although the detectors are enclosed in vacuum systems 
and isolated from vibrational, acoustic and electromag- 
netic disturbances, their typical output data contains a 
larger number of transient noise events (glitches) with 
higher amplitude than expected from Gaussian processes 
alone. Each observatory is equipped with a system of 
environmental and instrumental monitors that are sen- 
sitive to glitch sources but have a negligible sensitiv- 
ity to gravitational waves. These sensors were used to 
identify times when the detector output was potentially 
corrupted piH^ . We grouped these times into two 
categories: periods with well-understood couplings be- 
tween non-gravitational-wave sources and detector out- 
put, and periods when a statistical correlation was found 
but a coupling mechanism was not identified. In our 
primary search — which included the identification of 
gravitational-wave candidates and the calculation of up- 
per limits — we removed (vetoed) times that fell in ei- 
ther of the two categories from the analysis, along with 
any coincident events that occurred during these periods. 
We also performed a secondary search for possible loud 
candidate events, in which only the times with known 
couplings were vetoed. 

Approximately 10% of the data, designated play- 
ground, was used for tuning and data quality investiga- 
tions. This data was searched for gravitational waves, 
but not used in calculating upper limits. After all vetoes 
were applied and playground time excluded, there was 
0.09yr of HILIVI time, 0.17yr of HlLl time, O.lOyr of 
HlVl time, and 0.07 yr of LlVl time, giving a total anal- 
ysis time of 0.43 yr. 

A substantial change from the analysis procedure of 
[llj was that data were analyzed in two-week blocks 



with a latency of two to four weeks, to allow for feed- 
back of information to ongoing detector characterization 
efforts and to improve data quality. Thus, during the 
search many new vetoes were introduced resulting from 
improved understanding of the detectors. However, sig- 
nificant numbers of delta-function-like glitches with large 
amplitudes remained unvetoed in the LIGO detectors. 
These were found to cause artifacts in the matched filter 
output over a short time surrounding the glitch: thus, 
during the latter stages of the search, 8 s of time on ei- 
ther side of any matched filter SNR exceeding 250 was 
vetoed. Times removed from the primary search by this 
veto were still examined for possible loud events. 



IV. BLIND INJECTION RECOVERY 

The search pipeline described above identified a 
gravitational-wave candidate occurring on 16 September 
2010 at 06:42:23 UTC, with pc ^ 12.5 in coincidence be- 
tween the two LIGO detectors in the middle mass bin 
3.48 < M/Mq < 7 AO. The highest matched-filter SNR 
obtained in the search was 15 at = 4.7 in HI 
and 10 at = 4.4 M© in LI. This difference in SNRs 
is consistent with typical differences in antenna response 
factors for these differently-oriented detectors. Virgo was 
also operating at the time of the event, but its sensitiv- 
ity was a factor of approximately four lower than the 
LIGO detectors; the absence of a signal in Virgo above 
the single-detector SNR threshold of 5.5 was consistent 
with this fact. In the LIGO detectors, the signal was 
louder than all time-shifted HlLl coincident events in 
the same mass bin throughout S6. However, with only 
100 time-shifts, we could only bound the FAR to < 1/23 
years, even when folding in all data from the entire anal- 
ysis. To obtain a better estimate of the event's FAR we 
performed all possible multiples of 5-second time-shifts 
on four calendar months of data around the event, corre- 
sponding to an effective analysis time of 2.0 x 10^ years. 
We found five events with a value of pc equal to or larger 
than the candidate's, as shown in Figure (3] These five 
events were all coincidences between the candidate's sig- 
nal in HI and time-shifted transient noise in LI. When 
we excluded 8 seconds from around the event's time in the 
background estimation, we found no background events 
with Pc greater than the candidate and we obtained a sig- 
nificantly different background distribution, also shown 
in Figure [3j 

Including the events at the time of the candidate in 
the background estimate, the FAR of the event in the 
3.48 < M/Mq < 7.40 mass bin, coincident in the LIGO 
detectors, was estimated to be 1 in 4 x 10'* years. Since 
this event occurred in HILIVI time during VSR3, only 
two event types were considered: HlLl double- coincident 
events and HILIVI triple-coincident events. This re- 
sulted in a trials factor of 6 (accounting for the three 
mass bins and two coincidence types) and a combined 
FAR of 1 in 7, 000 years. The false alarm probability 
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FIG. 3: The cumulative rate of events with chirp mass 
3.48 < M/Mq < 7.40 coincident in the HI and LI detec- 
tors, seen in four months of data around the 16 September 
candidate, as a function of the threshold ranking statistic p^. 
The blue triangles show coincident events. Black dots show 
the background estimated from 100 time-shifts. Black crosses 
show the extended background estimation from all possible 
5-second shifts on this data restricted, for computational rea- 
sons, to only the tail of loudest events. The gray dots and 
crosses show the corresponding background estimates when 8 
seconds of data around the time of the candidate are excluded. 
Gray shaded contours show the 1 — 5a (dark to light) con- 
sistency of coincident events with the estimated background 
including the extended background estimate, for the events 
and analysis time shown, including the candidate time. This 
event was later revealed to have been a blind injection. 



of this event in this analysis, over the 0.47 yr of coinci- 
dent time remaining after all vetoes were applied, was 
7 X lO"'^. 

The detectors' environmental monitoring channels 
record data from seismometers, accelerometers, micro- 
phones, magnetometers, radio receivers, weather sensors, 
and a cosmic ray detector. Injections of environmen- 
tal signals and other tests indicate that these channels 
are much more sensitive to environmental signals than 
the gravitational wave readout channels are. Arrays of 
these detectors were operating and providing full cover- 
age at the time of the event, and did not record envi- 
ronmental signals that could account for the event. En- 
vironmental signal levels at our observatories and at ex- 
ternal electromagnetic weather observatories were typi- 
cal of quiet times. Mechanisms that could cause coinci- 
dent signals among widely separated detectors — such 
as earthquakes, microseismic noise due to large weather 
systems, and electromagnetic disturbances in the iono- 
sphere [231211 — were therefore ruled out. 

A loud transient occurred in LI 9 seconds before the 
coalescence time of the signal. That transient belonged to 
a known family of sharp (~ 10 ms) and loud (SNR 200- 
80000) glitches that appear 10-30 times per day in the 



output optical sensing system of this detector. Since the 
candidate signal swept through the sensitive band of the 
detector, from 40 Hz to coalescence, in less than 4 sec- 
onds, it did not overlap the loud transient. Studies, in- 
cluding re-analysis of the data with the glitch removed, 
indicated that the signal was not related to the earlier 
instrumental glitch. No evidence was found that the ob- 
served signal was associated with, or corrupted by, any 
instrumental effect. 

Following the completion of this analysis, the event 
was revealed to be a blind injection. While the analy- 
sis groups did not know the event was an injection prior 
to its unblinding, they did know that one or more blind 
injections may be performed during the analysis period. 
Such blind injections have been carried out before: see [1] 
for the results of a blind injection performed in a previ- 
ous run. This event was the only coherent CBC blind 
injection performed during S6 and VSR2 and 3. The in- 
jection was identified as a gravitational-wave candidate 
with high probability, and the blind injection challenge 
was considered to be successful [3]. 

In order to more accurately determine the parameters 
of the event prior to the unblinding, we performed coher- 
ent Bayesian analyses of the data using models of both 
spinning and non-spinning compact binary objects |27l - 
[5T] . These analyses showed evidence for the presence of 
a weak signal in Virgo, consistent with the signal seen 
by the two LIGO detectors. The strength of a signal 
in Virgo is an important input to the localization of a 
source in the sky. Parameter estimates varied signifi- 
cantly depending on the exact model used for the grav- 
itational waveform, particularly when we included spin 
effects. However, conservative unions of the confidence 
intervals from the different waveform models were con- 
sistent with most injected parameters, including chirp 
mass, time of coalescence, and sky location. In addi- 
tion, the signal was correctly identified as having at least 
one highly-spinning component with the spin misaligned 
with the angular orbital momentum. We will describe the 
details of parameter estimation on this and other CBC 
injections in a future paper (in preparation). 



V. SEARCH RESULTS 

After the event was revealed to be a blind injection the 
data containing it was removed from the analysis. With 
the injection excluded, there were no gravitational-wave 
candidates observed in the data. Indeed the search re- 
sult was consistent with the background estimated from 
time-shifting the data. The most significant event was an 
LlVl coincidence in LlVl time with a combined FAR of 
1.2 yr~^. The second and third most significant events 
had combined FARs of 2.2 yr^^ and 5.6 yr^^, respec- 
tively. All of these events were consistent with back- 
ground: having analyzed ~ 0.5 yr of data, we would ex- 
pect the loudest event to have a FAR of 2 ± 2yr~^. Al- 
though no detection candidates were found, a detailed 



9 



investigation of the loudest events in each analysis period 
was performed, to improve our understanding of instru- 
mental data quality. 



VI. BINARY COALESCENCE RATE LIMITS 

Given the absence of gravitational-wave signals, we 
used our observations to set upper limits on coalescence 
rates of BNS, BBH, and NSBH systems. We used the 
procedure described in [32H34] to compute Bayesian 90% 
confidence level upper limits on the coalescence rate 
for the various systems, making use of previous results 
[H Uni m] as prior information on the rates. 

The rate of binary coalescences in a spiral galaxy is 
expected to be proportional to the star formation rate, 
and hence blue light luminosity, of the galaxy [35. . Previ- 
ous searches [H UHl |TT] presented upper limits in terms of 
blue light luminosity, using units of L^g'^yr"^, where one 
Lio is 10^" times the solar blue light luminosity. There 
are, however, numerous challenges to evaluating the up- 
per limit as a function of luminosity, not least due to 
the large uncertainties in both the luminosity of and dis- 
tance to nearby galaxies, as well as the lack of a com- 
plete galaxy catalogue at larger distances [311 ES] ■ On 
large scales (greater than ^ 20 Mpc) , the luminosity per 
unit volume is approximately constant; consequently the 
analysis can be simplified by reporting upper limits per 
unit volume per unit time. During the current analysis, 
the sensitivity of the detectors to the systems of interest 
(as shown in Figure [5]) was sufficiently large that we could 
assume signals were uniformly distributed in volume. We 
therefore quote upper limits in units of Mpc^'^yr^^. To 
incorporate the previous results as prior distributions, we 
converted from Lio to Mpc'^ using a conversion factor of 
0.02 Lio per Mpc^ [35 . 

We estimate the volume to which the search is sensi- 
tive by reanalyzing the data with the addition of a large 
number of simulated signals ("software injections") in or- 
der to model the source population. Our ability to de- 
tect a signal depends upon the parameters of the source, 
including the component masses, the distance to the bi- 
nary, its sky location, and its orientation with respect 
to the detectors. Numerous signals with randomly cho- 
sen parameters were therefore injected into the data. To 
compute the sensitive volume for a given binary mass, we 
perform a Monte Carlo integration over the other param- 
eters to obtain the efficiency of the search — determined 
by the fraction of simulated signals found louder than 
the loudest foreground event — as a function of distance. 
Integrating the efficiency as a function of distance then 
gives the sensitive volume. 

We consider several systematic uncertainties that limit 
the accuracy of the measured search volume and therefore 
the upper limits [10] : detector calibration errors (conser- 
vatively estimated to be 14% in sensitive distance com- 
bined over all three detectors and over the entire ob- 
servational period, and a 2% bias correction), waveform 
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TABLE I: Rate upper limits of BNS, BBH and NSBH coales- 
cence, assuming canonical mass distributions. Dhorizon is the 
horizon distance averaged over the time of the search. The 
sensitive distance averaged over all sky locations and binary 
orientations is Davg — -Dhorizon/2.26 |36) . The first set of up- 
per limits are those obtained for binaries with non-spinning 
components. The second set of upper limits are produced 
using black holes with a spin uniformly distributed between 
zero and the maximal value of Gm? /c. 



errors (taken to be a one-sided 10% [32] bias towards 
lower sensitive distance) , and Monte Carlo statistical er- 
rors (3-5% in sensitive volume) . We convert the sensitive 
distance uncertainties to volume uncertainties, and then 
marginalize over the uncertainty in volume to obtain an 
upper limit which takes into account these systematic 
uncertainties |32| . 

In Table [ij we present the marginalized upper limits at 
the 90% confidence level assuming canonical mass dis- 
tributions for non-spinning BNS {mi — m2 — 1.35 ± 
0.04Mq), BBH (mi = ms = 5 ± IM©), and NSBH 
(toi = 1.35 ± 0.04 Mq, TO2 = 5±1M0) systems. We also 
compute upper limits as a function of total mass M, us- 
ing an injection population distributed uniformly over M 
and uniformly over mi for a given M. For NSBH systems 
we present the upper limit as a function of black hole 
mass, keeping the neutron-star mass fixed in the range 
1 — 3M0. These are presented in Figure ]4] Figure ]5] com- 
pares the upper limits obtained in this analysis (dark gray 
regions) to limits obtained in our previous searches up to 
S5/VSR1 [4 (light gray region) and to astrophysically- 
predicted rates (blue regions) for BNS, NSBH, and BBH 
systems. The improvement over the previous limits is up 
to a factor of 1.4, depending on binary mass; this reflects 
the additional observation time and improved sensitivity 
of the S6/VSR2/VSR3 data with respect to all previous 
observations. 

Although we searched with a bank of non-spinning 
templates, we compute upper limits for NSBH and BBH 
systems in which one or both of the component masses 
are spinning. These results are also presented in Ta- 
ble |T] We did not compute upper limits for spinning 
BNS systems because astrophysical observations indicate 
that neutron stars cannot have large enough spin to sig- 
nificantly affect waveforms observable in the LIGO fre- 
quency band [371 ISH]- Black hole spins were uniformly 
distributed in both orientation and magnitude, S, with 
S constrained to the range < S < Grn^/c, and m is 
the mass of the black hole. As can be seen in Table [T| 
the spinning upper limits are ~ 16% larger than non- 
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2.0 5.0 8.0 11.11 11.0 17.0 20.0 25.0 BNS NSBH BBH 

Total Mass (Mq) 




3.0 8-0 13-0 18.0 23.0 

Component Mass (M0) 



FIG. 4: The marginalized upper limits as a function of mass. 
The top plot shows the limit as a function of total mass M, 
using a distribution uniform in mi for a given M. The lower 
plot shows the limit as a function of the black hole mass, with 
the neutron star mass restricted to the range 1 — SM©. The 
light bars indicate upper limits from previous searches. The 
dark bars indicate the combined upper limits including the 
results of this search. 



FIG. 5: Comparison of CBC upper limit rates for BNS, NSBH 
and BBH systems. The light gray regions display the upper 
limits obtained in the S5-VSR1 analysis; dark gray regions 
show the upper limits obtained in this analysis, using the S5- 
VSRl limits as priors. The new limits are up to a factor of 
1.4 improvement over the previous results. The lower (blue) 
regions show the spread in the astrophysically predicted rates, 
with the dashed-black lines showing the "realistic" estimates 
[5]. Note: In |5], NSBH and BBH rates were quoted using a 
black-hole mass of 10 M©. We have therefore rescaled the S5 
and S6 NSBH and BBH upper limits in this plot by a factor 
of (Ads/Adio)^^^, where Alio is the chirp mass of a binary in 
which the black hole mass is 10 M© and AI5 is the chirp mass 
of a binary in which the black hole mass is 5 M© . 

3 X 1G~* Mpc~^yr~^ with at least an order of magnitude 
uncertainty in either direction [S]. In all cases, the upper 
limits derived here are two to three orders of magnitude 
above the "realistic" estimated rates, and about a fac- 
tor of ten above the most optimistic predictions. These 
results are summarized in Figure [5] 

VII. DISCUSSION 



spinning. Signals from spinning systems are recovered 
with a worse match to our templates since we use a non- 
spinning template bank. 

While the rates presented here represent an improve- 
ment over the previously published results from ear- 
lier LIGO and Virgo science runs, they are still above 
the astrophysically predicted rates of binary coalescence. 
There are numerous uncertainties involved in estimat- 
ing astrophysical rates, including limited numbers of 
observations and unknown model parameters; conse- 
quently the rate estimates are rather uncertain. For 
BNS systems the estimated rates vary between 1 x 10~* 
and 1 X 10~^ Mpc~'^yr~^, with a "realistic" estimate 
of 1 X 10-*^Mpc"^yr-i. For BBH and NSBH, realis- 
tic estimates of the rate are 5 x 10~^ Mpc^'^yr"^ and 



We performed a search for gravitational waves from 
compact binary coalescences with total mass between 2 
and 25 M0 with the LIGO and Virgo detectors using 
data taken between July 7, 2009 and October 20, 2010. 
No gravitational waves candidates were detected, and we 
placed new upper limits on CBC rates. These new limits 
are up to a factor of 1.4 improvement over those achieved 
using previous LIGO and Virgo observational runs up to 
S5/VSR1 [4], but remain two to three orders of magni- 
tude above the astrophysically predicted rates. 

The installation of the advanced LIGO and Virgo de- 
tectors has begun. When operational, these detectors will 
provide a factor of ten increase in sensitivity over the ini- 
tial detectors, providing a factor of ~ 1000 increase in 
the sensitive volume. At that time, we expect to observe 
tens of binary coalescences per year [S] . 

In order to detect this population of gravitational wave 
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signals, we will have to be able to confidently discrimi- 
nate it from backgrounds caused by both stationary and 
transient detector noise. It is customary [S] to assume 
that a signal with SNR 8 in each detector would stand 
far enough above background that we would consider it to 
be a detection candidate. The blind injection had some- 
what larger SNR than 8 in each detector, and we were 
able estimate a FAR of 1 in 7000 years for that event. 
Alternatively, consider a coincident signal with exactly 
SNR 8 in two detectors. Provided the signal is a good 
match to the template waveform (x^ ~ 1 in equation [l]) 
this corresponds to pc = 11.3. As can be seen from the 
extended background events with the blind injection re- 
moved in Figure [s] (light-gray crosses) , this gives a FAR 
of ~ 1 in 2 X lO'^ years in a single trial, or 1 in 3000 
years over all trials. Achieving similar-or-better back- 
ground distributions in Advanced LIGO and Virgo will 
require detailed data quality studies of the detectors and 
feedback from the CBC searches, along with well-tuned 
signal-based vetoes. We have continued to develop the 
pipeline with these goals in mind. For this analysis we 
significantly decreased the latency between taking data 
and producing results, which allowed data quality vetoes 
to be finely tuned for the CBC search. These successes, 
along with the successful recovery of the blind injection, 
give us confidence that we will be able to detect gravi- 
tational waves from CBCs at the expected rates in Ad- 
vanced LIGO and Virgo. 
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